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I. INTRODUCTION 

The one-dimensional(lD) quantum spin chains play a very important role in the strongly correlated systems 
and low-dimensional quantum magnetism and show many interesting behaviors. For example, the ground state 
of the spin-1/2 Heisenberg chain with antiferromagnetic couplings is expected to be a spin liquid rather than Neel 
ordered state due to the strong quantum fluctuations. The elementary excitations of such a system are usually 
described by the spinons which carry spin-1/2 rather than spin waves. In another hand, Haldane conjectured that 
the spin systems with half-integer spins have gapless excitation spectra, while those with integer spins have gapped 
spectra The spin-1/2 Heisenberg model was exactly solved by Bethe Q. By using Bethe's hypothesis, Yang 
and Yang solved the XX Z Heisenberg chain successfully @. Subsequently, Takhtajian and Faddeev developed 
the algebraic Bethe ansatz method and several spin chain models have been exactly solved [H [f| @, 0, 

Recently, much attention has been focused on the high spin systems because not only many peculiar quantum 
orders and exotic collective excitations can appear in this kind of systems, but also some materials such as 
CsVBr 3 , CsVCbj [l||, CSVI3 and AgCrP 2 Se [13] in nature can be modeled by the spin-3/2 chain quite well. 
With the developments of experimental technique of laser cooling and magnetic traps, atoms with high nuclear 
spin such as 87 Rb and 23 Na with spin-1, 132 Cs, 9 Be, 135 Ba, 137 Ba, and 53 Cr with spin-3/2 can be trapped. Using 
the Feshbach resonance techniques, one can tune the scattering lengths among the atoms, which make it possible 
to simulate the traditional solid systems with various interactions. These progress provide us a ideal plat form to 
study the physics in high spin systems and many interesting quantum phenomena are found [la , 1 19L [2fJ . 

It is well-known that the spin-s chain with the SU(2s + l)-symmetry can be solved exactly, for the Hamiltonian 
can be mapped onto the summation of permutation operators. Takhatajan and Babudjian found that besides 
the SU(3) integrable point, the spin-1 system can still be solved at a special S'J7(2)-invariant point [1, [|| . They 
also showed that the elementary excitation of the system is gapless. This motivate us to seek other integrable 
points in the high spin systems. In this paper, we show that besides the SU(4) integrable point, there is another 
exactly solvable model of the spin-3/2 chain with SO (5) symmetry. By using the nested algebraic Bethe ansatz 
method, we obtain the exact solutions of the model. Based on the exact solutions, several exotic excitations such 
as the heavy spinon with fractional spin 3/2, the neutral spinon with spin zero and the dressed spinon with spin 
1/2 are found, which are quite different from those of the ST/ (4) integrable spin chain. 

The paper is organized as follows. We introduce the model and its symmetry in Sec. [ITJ The nested algebraic 
Bethe ansatz approach for the model is shown in Sec. IIIII The thermodynamic properties of the system are 
analyzed in Sec. HVi The ferromagnetic and antiferromagnetic ground states are discussed in Sec. [V] The 
elementary excitations are given in Sec. IVTl Sec. I VIII is a brief summary. 



II. THE MODEL 



As mentioned above, two integrable models of ID spin-3/2 Heisenberg chains, i.e., the SL/(2)-invariant one and 
the St/ (4)-invariant one, have been found and solved. In this paper, we introduce another integrable spin-3/2 
chain model with SO(5)-symmetry. Our model Hamiltonian reads 



iY 



H = J J2 



— S, • Si+i - - (^Si ■ Si+i ; 



2 / 3 - 
- ySi ■ Si+i 



(l) 



where J is a coupling constant; Si is the spin-3/2 operator at site i, i = 1, 2, • • • ,N and N is the length of the 
system. Here, we adopt the periodic boundary condition, i.e., Sn+i = S±. 
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To show the SO(5) symmetry of our model, we introduce the following Dirac matrices 

r 1 = —-^=(s x Sy + SyS x ), r 2 = -^=(s z s x + s x s z ), r 3 = — j=(s z Sy + s y s z ) 7 

f 4 = s 2 -5/4, r 5 = i=(s 2 -s 2 ). (2) 

The 10 generators of the 5*0(5) Lie algebra can be expressed by the Dirac matrices as F a ' b = — |[r a ,F 6 ]. The 
explicit from of these generators are [21| 

F l, (2.3,4) _ ( a A 1,5 _ / -I 
~ \a )> 1 ~ \ I 

r 3,4 ;4)2;2l 3 = (9 0\ r2 , ( 3, 4l5) = / -ia j 



v o <r y ' yia o 

After some algebra, we find that the Hamiltonian (1) commutes with the generators Q, [H,T a ' b ] = 0. Thus the 
model has the SO(5) symmetry. Different from the ST/(4) integrable spin chain, there only exist three conserved 
quantities and the number of spins with individual components is no longer conserved. After some detailed 
analysis, we find that the following quantities are conserved 

Ji = N 3/2 + N 1/2 + N_ 1/2 + N_ 3/2 , 

■h=N V2 ~N_ m , (4) 
J3 = N 1/2 -N_ 1/2 . 

The i?-matrix of this model reads 

RabW = -^Yi P °» + p « b - + p *»> (5) 

where A is the spectral parameter; P^ b is the projection operator in the total spin-s channel and acts on the two 
coupled spin space V a (£> V b - Just as the ST/ (4) one, the non-trivial scattering processes only exist in the total 
spin-0 and 2 channels, but here the scattering strengths are different in these two channels. 
This i?-matrix satisfies the Yang-Baxter equation [1, 

R ab {X)R bc {X + ^RabW = Rbc(n)Rab(\ + v)RbcW- (6) 

In the frame work of quantum inverse scattering method(QISM), the Lax operators of the system are 

L n (X) = RonPon = 2A - 3i P ° b + P ° b + 2A - i P ^ b + P ^ b ' ^ 
where Vq is an auxiliary space and V n is the quantum space. Pq 71 is the permutation operator, which can be 
expressed by the projection operators as Pq u = — P(f n + Pn n — Pin + Pon- The monodromy matrix T is constructed 
by the Lax operators as 

T(A) = L N (\)L N _i(X) ■ ■ ■ Li(X). (8) 

The monodromy matrix T satisfies the Yang-Baxter relation 

R(X - /Lt)T(A) ® T(/i) = T(n) ® T(\)R(X - n), (9) 

Taking trace in the auxiliary space of T, we obtain the transfer matrix 

t(X) = tr T(A). (10) 

From the Yang-Baxter relation, one can prove that the transfer matrices with different spectral parameters 
commutate with each other, 

[t(A), *(/*)]= 0. (11) 

These transfer matrices are the infinite conserved quantities of this system. Thus the system is integrable. 
Taking the derivative of the logarithm of the transfer matrix, we arrive at the Hamiltonian (1) (23l. [24l| 



9J SlnT(A) 

H = — 



4 dX 



9 4jN = -3 J (P° i+1 + 3^ l+1 ) - ? JN, (12) 

A=0 



Here, we have used the fact that the Lax operator ([7]) degenerates into the permutation operator when the spectral 
parameter is zero, L ab (X)\x=o = Pab, and put a constant J into the Hamiltonian. The eigenvalue problem of 
the Hamiltonian is therefore turned into the diagonalization of transfer matrices. Suppose the eigenvalues of the 
Hamiltonian and the transfer matrix are E and A, then the eigenvalue E can be determined by the eigenvalue 
A as 



E = J_jdhiA{X) 



A r)\ 



qq 

i;N. (13) 
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III. NESTED ALGEBRAIC BETHE ANSATZ 



In this section we show that the system (1) can be solved exactly by using the nested algebraic Bethe ansatz 
method similar to that used for the spl (2\1) supersymmetric one [25j . 

The local vacuum state of the zth site is chosen as |0) ^ = |3/2) i = (1,0,0,0)*, where t means the transport. 
The Lax operator on the site i can be written into following matric from in the auxiliary space 



Li(\) 



( 



\ 



A(X) *-Bi(A) 

C\X) W\(X) 

C 2 {X) l D 2 (X) 

G(A) 'CffA) 



l 5 2 (A) l F(X) \ 

l Dl(X) 'B^iX) 

4 D|(A) l B* 2 {X) 

l C* 2 {X) i V(X) ) 



(14) 



where M, l D 3 k , l Bj, '£?*, l C, l Cj, l C*, l G : l V and l F are operators in the quantum space Vi. The Lax operator 
Li acting on the local vacuum state gives 



Li(\)\0) 



(l l Bi(A)|0) i B 2 (X) |0) t F{X)\0)\ 

/(A) |0> i S* 1 (A) |0> 

/(A) |0> *S* a (A)|0) 

V p(A) |0) J 



(15) 



where /(A) = 2A/(2A + i). The operators 4 Bi(A) and l B 2 (X) are 



%(A) 



/ 

m(A) 



\ -s(A) 



l i? 2 (A) 






m(A) 




s(A) 




(16) 



where m(A) = i/(2A + i) and s(A) = — 2A«/(2A + 3i)(2A + i). Thus the actions of % B\ and *i?2 get a flip of 
one and two spin quanta respectively. Obviously, the whole state space can be obtained by these two operators. 
This is quite different from the 5'C/(4) one in which | — 3/2)^ can not be reached with the flips made by % B\ and 
l B 2 . So that the construction of the eigenstates might only need these two operators in this model. This makes 
that only once nesting is needed in the algebraic Bethe ansatz of the model (1). In fact l F is also needed in the 
construction of the eigenstates but plays an assistant role and doesn't affect the number of nestings. 

The global vacuum state of the system (fTJ) is the production of local vacuum states |0) = 10^® |0} 2 <8> • ■ -® |0)jy. 
The monodromy matrix T acting on the global vacuum state gives 



T N {X) |0) - L N {X)L N ^{X) ■ ■ ■ Lx(X) |0) 

/I, E^ffiPO 4 Si(A), E^|(A) *B 2 (A), 



o, 


f N W, 


o, 


o, 


o, 




\°. 


0, 





E i5l 3 (A) ^(A) 
E t 9tW l B* 2 (X) 
V 



|o>, 



(17) 



where g\ are some coefficients, and * represents the nonzero element. The matric form of the monodromy matrix 
T in the auxiliary space is 



T(A) 



/ A{X) Bx{X) B 2 (X) F(X) 

CHX) D\(X) D\{X) B*\X) 

C 2 (X) D\{X) D 2 {X) B* 2 {X) 

\G(X) C*(X) C* 2 {X) V{X) 



whereA, B h B*i , G, C\ C* 



D J 



F and V are operators in the Hilbert space V\ <£> V 2 



(18) 



Vn- Acting these 



elements on the vacuum state, we obtain 



a^(A)|0) = 0, C;(A)|0) = 0, G(A)|0)=0, 
A(X) |0) = |0) , Dl(X) |0) = Dl{X) |0) = f N (X) 



D\{X)\Q) = 
3), V(X) 



D?(A)|0)=0, 
))=P N (X)\Q). 



(19) 
(20) 



The operators Cj, C* , G, D\ and D 2 acting on the vacuum state |0) gives zero and the operators A, D\, D 2 
and V acting on the vacuum state |0) give the eigenvalues. The elements Bj, F and B* J acting on the vacuum 
state |0) will generate other spin-flipped states, and thus can be regarded as the generating operators of the 
multi-particle eigenstates of the transfer matrix. 

Now we turn to the eigenvalue problem of transfer matrix T. From the definition p8[) . the transfer matrix 
can be expressed in the form of 
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The eigenvalues of T can be determined by the eigenvalues of operators B, D{, D\ and V. These operators are 
shown to be diagonalized in vacuum state |0) in Eq. ([20]) . Obviously, 1-00) = |0) is an eigenstate of t(X), and the 
corresponding eigenvalue is 

A (X)=[l + 2f N (X)+p N (X)]. (22) 

Here \ipo) 1S called zero-particle state. To construct the other eigenstates of the transfer matrix, the generating 
operators Bj, B*^ and F should be used. Bj and F are enough to generate these states as shown bellow. Assume 
that the eigenstates of the transfer matrix have the form \ip) = %p(Bj,F) |0). To go further, we need need the 
commutation relations of A{X)B a {p)^D b {X)B c (^) , V(X)B a (fx), B(X)F(fj,), A£(A)F(/i) and V(X)F(/j,). From the 
Yang-Baxter relation © we obtain |25j 

Til 

A(X)B a (fj.) = -jB a (X)A(n) + jB a (fi)A(X), (23) 
D a b (X)B c (p) = jB e (n)D a d (X)r(X - fx)f c - jB b {X)D a c {») 



—B* a A([i) + -F(X)C a {fi) - -F(/x)C°(A) 
P P P 



6c (24) 



Where p(X) = 4A(A + i)/(2X + 3i)(2A + i), £ is a vector of 2 <£> 2-dimension, £ = (0, 1, —1, 0), r is the matrix 



r(A) = 



1 0\ 

6(A) o(A) 

o(A) b(X) 

1 / 



(25) 



where a(X) = A/ (X+i), 6(A) = i/(X+i). Here the parameter (X — fi) is omitted for short and the parameter (fi— A) 
is also omitted by adding a tilde above the corresponding function. The commutation relation of V(X)B a (/i) is 



(26) 



K(A)B„(aO = f -BMV{\) + -F{n)C* a {X) - -F(X)C*M - -^ bc B* b (X)D c M, 
P P P P 

where n(A) = ic{AX + 3i)/(2A + 3i)(2A + i). 

A{X)F{^) = If(h)A(\) - "F(\)A(n) + -MX)B b {^)C\ (27) 
P P P 

D a b (X)F(») = (1 - ?£)FQi)Dt(\) + ^F(X)D%(P) 

+ jiB* a (X)B b (ri-B b (X)B* a (n)}. (28) 

The commutation relations of D"F are contained in them. 

V(X)F(n) = -F(n)V(X) - -F(X)V(fi) ~ -UB* a (X)B* b (ri, (29) 
P P P 

For C a and C* appear in the left hand of (|26p. the commutation relation of them with the operators B a and F are 
also needed in the discussion. For only one and two-particle cases are discussed in detail here, the commutation 
rules used in this article are 

77 

C* a {X)B b {n) = -B a ^)C* b {X) + B b (f,)C* a (X) 

+ -[A^)V(X) - F(n)G(X)]U - -B a (X)C* b (p) - -UD c a (X)Dt(fi), (30) 
P P P 

777 

C a (X)B b {^) = B b (n)C a (X) + j[Afa)DS(\) - A(X)D a b (p)}, (31) 

B a (X)B b (f,) = B c (fi)B d (X)r(X - + -[F(\)A(jj>) - F(^)A(X)]U, (32) 

P 

B a (X)B*\fi) = B* b (fx)B a (X) + ZiF^D'jX) - F(X)D b M]U- (33) 
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Additionally, in the discussion of the three-particle wave functions, the following commutation relations are 
needed, 



B a (X)F(») = ■jF(ji)B a (X) - jF(X)B a (»). 



(34) 



From the commutation rules listed above, the eigenstates and eigenvalues of the transfer matrices can be 
discussed. The one particle state can be defined by Bi operator acting on the vacuum state as usual, 



\ip 1 (X 1 ))=B a (X 1 )W a \0), 
where W a the coefficients of B a (Xi). The transfer matrix acting on the one particle state arrives 
t(X) |^(Ai)) = [A(X) + D n (X) + D 22 (X) + T/(A)]B a (A 1 )^ Q |0) 

1 ■ f{X - Xl \ N (X) + J^4 T [l + a(X-X 1 )}\\M^)) 



(35) 



/(Ai - A) p(X - Ai) v ' /(A - Ai) 

f N (Xx)\B a (X)W a \0) 



m(Ai — A) m(X — Ai) N 



/(Ai - A) /(A -AO 



s(A-Ar) 



P(A-Ai) 



/at(Ai) - 1 CafcS* 6 (A)^ a |0) 



(36) 



If the one particle state is a eigenstate of the transfer matrix, terms including B a {X)W a |0) and B£(X)£ a bW a |0) 
should be canceled with each other. The condition that the unwanted terms cancel with each other gives 



/ L (Ax)-l = 0. 



(37) 



The Eq. ([37)) is the Bcthe ansatz equation. If the parameter Ai satisfies the Bethe ansatz equation (|37l) . the one 
particle state (j35jl is an eigenstate of the system. The eigenvalue A\(X, Xi) of the transfer matrix at one-particle 
state is 



Ai(A,Ar) 



1 



/(A ~ A V(A) 



/(Ai-A) p(A-Ai) 



/(A -AO 



[1 + a(X - AO]- 



(38) 



To generate the two particle states, operator F is needed additionally. The two-particle eigenstate is assumed 
to be 



|^ 2 (Ai, A 2 )) = B ai (A0-B a2 (A 2 )I^ a2ai |0) + ^(Ax, A 2 )F(A0Ca 2Ql ^ a2ai |0) 



(39) 



where h(Xi,X 2 ) is a undetermined function. Acting the transfer matrix ^at(A) on the assumed state (|39p . we 
obtain 



MA)h&2(Ai,A 2 )) = [A(X) + D{(\) + D 2 2 (\) + V(\)] |^(Ai,A 2 )) 

= + + |vl) + |V|) + |V> 2 ) + IV'I) + k 6 2 ) + 



(40) 



Here, \^ 2 ) denotes the eigenstate which including the operators B ai (Xi)B a2 (A 2 ) and F (Ai). V 1 !)) l^l); 

IV'I); (V'i)) l^f) and (V'l) denote the unwanted terms including B* i (X)B* 2 (Xi), B m (X)B* n (Xi), B ai (X)B a2 (X 2 ), 
B^ n (X)B a2 (X 2 ) 1 B ai (X)B a2 (Xi), B* m (X)B c (\i) and F(X) respectively. The unwanted terms should be canceled 
with each other, which gives the form of the function h and the Bethe ansatz equations. The unwanted terms 
l^!) and \ip 2 ) can be explicitly expressed as 



a(A-Ai) 
p(A-A0 

m(A — AO 



/(A -AO 



h(X!,X 2 )- 
h(X 1 ,X 2 )- 



fl(Ai-Aa) 



p(Xi - A 2 ) 
s(Ai - A a ) 



If the function h satisfies 



h{X u X 



p(Xi - A 2 ) 

«(Ai - A 2 ) 



^ lb2 W blb ^ aia2 B* a ^X)B* a ^X 2 ) |0) , 
Ua 2 W a ^B m (X)B: n (X 2 )\0). 



p(Xi - A 2 ) ' 



(41) 
(42) 

(43) 



G 



the terms Jf/^) an d (V'i) are canceled with each other. This constraint determines the values of h introduced in 
the two-particle eigenstate ([M]). The unwanted terms j^f) an d are 



77l(Al — A) 



/(Ai - A)/(A 2 - Ai) 
mCA-AO/^CAx) 



"2(11 U/i>lf>2 



B Ol (A)B„ 2 (A 2 )|0> 



s(A-Ax) 



p(A - Ai)/(A 2 - Ai) 



B^(A)B Oa (A 2 )|0). 



(44) 



(45) 



After some tedious calculations, we find that if the parameter Ai and A 2 satisfy the following Bethe ansatz 
equations 



/ JV (A 1 )/(A 2 -A 1 ) 1 ^ nAl A2 Vfc 2 ^ ■> 

the terms i/>|) an d ^2) are a ^ so canceled with each other. The unwanted term J^f) is 



(46) 



i>l) = \h(X u X 2 ) 



£(Ai-A) coia2T/r/6l62 , m(Ai-A)m(A 2 -Ai) T;r/aiaa 



p(A x - A) 
m(A 2 - A) 



£h h £ aia2 W 0l ° 2 H — — — —W a 

fetlt + /(Ai - A)/(A 2 - Ax) W 



/(A - Ai)/(Ai - A 2 ) 
Wc find that if the parameter Ai and A 2 satisfy Eq 

1 /(A 2 - x i) w a ia 



/(A 2 - A)/(Ai - A) 



and 



/ N (A 2 )/(A!-A 2 ) 



(47) 



(48) 



The unwanted term Ji/jf) is zero. The Eqs. (|46|) and |48|) are the two-particle Bethe ansatz equations, which can 
be written into a uniform form 



}_ /(Ai - Aj) a 



r(Ai) /(A, - A,) 



r{*i-* i )%£W*'», = 1,2, 



(49) 



The function W aia2 can be determined by the nested algebraic Bethe ansatz. Together with Eq. (j43]) . the 
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following unwanted terms J'i/'f) an( i I^I) can be canceled, 

-s(A-A 1 )m(A 1 -A 2 ) /(A - A x ) s(A - A 2 ) 



1^) 



p(A-Ai) /(A1-A2) p(A-Ai)p(A- A 2 ) 



/ iV (A 2 )^ mQ2 iy c 



s(A - Ai) m(A 2 - Ai) 
'p(A-Ai) /(A2-A1) 

x w aia2 



s{X-X 2 ) 



/(A-Ax)p(A-A 2 ) 



S™(A)S C (A 1 ) |0) + h(X 1 ,X 2 )^-^-^ 2ai W a ^B* m (X)B m (X 1 ) |0) , 



/(A-Ai) 
fl(Ai-Aa) 



p(Ai - A 2 ) 



n(A-Ai) 
P(A-Ai) 

1 m(A 2 - A) s(Ai - A) .n(Ai-A) 

+ ft(Ai, A 2 J 



P JV (A 1 )F(A)6, 1 « a W° iaa |0> 



/(Ai - A) /(A 2 - A) P (X 1 - A) 



'p(Ax-A) 



p W aia2 

C,a 1 a 2 W 



xF(X) |0> 



n(X - A x ) s(Ai - A 2 ) m(A - Ai) s(A - A 2 ) 



p(A - Ai) p(Ai - A 2 ) p(X - Ai) p(A - A 2 ) 



xr i0! f(A) |0) + 2/i(Ai,A 2 ) 



m 2 (A - Ai) 
/ 2 (A-A X ) 



r(A!)^ 2ai ^ aia2 F(A)|0) 



+ 



m(X — Ai) s(A — A 2 ) 



/(A- 


Ai) p(X- 


-Ax) p(A- 


A 2 )^ 


s(A- 


Ai) m(A 


- Ai) m(Ai 


-A 2 ) 


/(A- 


Ai) p(A- 


-Ai) /(Ax 


-A 2 ) 


1 


m(A 


- A 2 ) s(Ai - 


-A) , 
-r 



Nbl&2 fN I 



u ai [r(Xi)-r(X2)]w a ^F(x)\o) 



KA - AO^ / JV (A 2 )C &lb2 ^^F(A) |0) 



/(A — Ai) p(A-A 2 ) p(Ai-A) 
With the above results, the explicit form of the two-particle eigenstate of the system is 

1 



1*2) 



/(A-A0/(A-A 2 ) N 

P ( A ) 



/(Ai - A)/(A 2 - A) p(Ai - A)p(A 2 - A) 
f N W 



+ 

+ { (Ai,A ; 
1 



■ A l )/(A-A 2 ) W A _ Xl >^ X X ^J BaMBaMW™ |0) 
1 p N (X) ( m 2 (A-Ai)^ 



p(Ai-A) p(A-Ax) 



1 - 



/ 2 (A-A X ) 



m(A 2 — A) s(Ai — A) ( s(A — A 2 ) m(A — Ai) N 



/(Ai - A) /(A a - A) p(Aj - A) p(A - A 2 ) /(A - Ai) 



2f N (X) 

p n W 



/(A — Ai)p(A — Ai) /(A-A 2 ) 
1 



/(A -Ax) 



m(A- A 2 ) s(Ai - A) q 2ain w 
/(A-A 2 )p(A 1 -A) r6 ^ (A " Al) W (A2) 



s(A-A 2 )m(Ai-A) 6 ^ 
p(A-A 2 ) p ( Al -A) Ai)^ 2( j 2 / (A) 



F(Ai)W a2Ql |0). 

Considering Eqs. (j43|) and |49|) . the eigen-equation of the transfer matrix becomes 



*(A)iMA l! A 2) )=n 7 ^ + p-(A)n^M + n. j 



l\f(Xi-X) * v '^piX-X,,) AA/(A-A,) 
xKA - Ai)«^r(A - A 2 )^^ ld2 [Sa 1 (A 1 )B a2 (A 2 ) + /i(Ai, A 2 )F(Ai)£ a2ai ] |0) 
The eigen- value of two-particle state is 

2 , 2 



A 2 (A,{A 1 ,A 2 }) = J 



L/(A,-A) P 1 Ul^A-A,) 



r(A)n 



L/CA^^^MAr.A,}), 



(50) 



(51) 



(52) 



(53) 



(54) 
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From the above discussions, we see that the construction of eigenstates of the 5'0(5)-invariant quantum spin 
chain is quit different from that of the iSE/(4)-invariant one. The spin-flipped operators are very complicated. The 
symmetry analysis of the states is helpful to construct the eigenstates [25j . Now, we construct the many-particle 
eigenstates. We define ipn as the creation operators of the n-particle state for convenience 

\i> n {{\i})} = ip n W \0) , (55) 

where W = {W X ,W 2 ,--- ,W n ) T are some vectors. Obviously, tp = l for \tp ) = |0). From Eqs. ([35]) and (|39j) . 
we know the creation operators of one and two-particle eigenstates are 

$i = B(\ 1 ), $ 2 (X 1 ,X 2 ) = B(X 1 )®B{X 2 ) + h{X u X 2 )F(X 1 )£ (56) 

From the commutation relations (f3"2"| . the relation of A 2 ) and 4>(X 2 , Ai) is 

^2(A2,Ai) = -0 2 (Ai, A 2 )r(A 2 - Ai), ip 2 (Xi,X 2 ) =tp 2 (X 2 ,Xi)r(X 1 - X 2 ). (57) 

After a more detailed analysis, we find that the three-particle states should have the same symmetry. Assume 
the three-particle eigenstates as 

MXi, X 2 , A 3 ) = S(Ai) <g> MX 2l A 3 ) - ft(Ai, X 2 )F{X!)t® ^(A 3 ) 



/(A 3 - A 2 ) 



which satisfies 



h(X 1 ,X 3 )F(X 1 )^ <g> Vi(A 2 ) 77 r ^^(Aa - A 3 ), (58) 

J{A 2 — A 3 j 



|V'3(A 2 , Ai, A 3 )) = |^' 3 (Ai, A 2 , A 3 )) r i2 (Ai - A 2 ), 

|V>3(Ai,A 3 ,A 2 )) = |^ 3 (A 1 ,A 2 ,A 3 ))r 23 (A 2 - A 3 ). (59) 



Borrowed the ideas in [25j , the generator of multi-particle state is assumed as 
$Mi(Xi, A 2) - ■ ■ ,A Ml ) = B(Xi) ® Vmi-i(A2, A 3l - ■ ■ ,A Ml ) 

Ml n 

- ^(Ai)!*® h(Xi - Xj) [] 



f(X k - A,) 

j= 2 k=1,k^ 3 J y J > 

x ipM 1 ~ 2 (X 2 , ■ • • , Aj_i, • ■ ■ , Aj+i, AmJ rjb,fe + i(Aft - Aj), (60) 

fe=2 

which satisfies 

^Mi(Ai, A 2 , • • • , X Ml ) = V'a/i(Ai, • • • , X j+ i,Xj, • • • , \ Ml )rj t j + i(Xj - X J+ i). (61) 
By using the same process as that for the two-particle case, we obtain the following Bethe ansatz equations 

Mi 



fjXj - Xj) 



— VU-ai-a Ml _ TT 

J I\A 11 



x >l&i-i(Ai, {A i+1) • • • X Ml , A l7 • • • A^x})^-^ , (62) 



Mi 



A Ml (\{*i}) = II tn 1 x/ S(A,{Ai», (63) 
fc=i ^ Afc " A ) 

where 7Lj^(A, {Ai}) is the eigenvalue of a series production of r matrices 

r(X - Ai)££r(A - A 2 )$° ■ ■ -r(A - X n )t : ^ d W d ^ = (A, {A,}) W"'"" 1 . (64) 

Now, we diagonalize the eigen-equation (f6"4l . One can easily check that the nested R matrix satisfies the 
Yang-Baxter equation ([6]). In fact, it is the 6-vertex i?-matrix in space <g> with a two dimensional 
space. The corresponding Lax operator is 
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The nested monodromy matrix is denned as 

zffi(A, {A,}) = L^ Mi (A - AOL^^^A - A a ) • • • lQ{X - X Ml ), 
which satisfies the Yang-Baxter relation ©. The nested transfer matrix is 

4)(A) = tr r£(A). 

The transfer matrices with different spectral parameters commutate with each other 

[i«(A),4>)]=0. 



(66) 
(67) 
(68) 



Using the standard Bcthc ansatz method for the six-vertex model, we obtain the eigenvalue of nested transfer 



matrix t$ as 



M, 



M-2 



<(A,{^})=ri' 



Mi 

n«(A-Ai)n 

i=l 



where the parameter /ij should satisfy the following Bethc ansatz equation 

ii ^4=n«<«-A,). 



Substituting Eq. 



t . o(/x !; - Aj) 

4=1, X ' 1=1 

into (f6"2")l and we obtain the eigenvalues of T/v(A) as 

Mi 1 Mi 



(69) 



(70) 



f 



A M iM 2 (A,{Ai},{^}) = n /( 



p n WI 



Ml 



i 



Mi 



M 2 



/(A - Aj) 
Li p(A - A t ) 

M 2 



The Bethe ansatz equations are Eq. ([70]) and 

1 



]a(A-A fc )J 



1 



i—k 



\a(\-Hj) ^a(Mi-A) 



n 



I^Mi) 



/ L (Aj) 



rf M - Aj) ft 
11 f(\_. _ 11 



/(Aj - A ; ) — Aj)' 



(71) 



(72) 



Putting A; — > Aj — i/4, ^ — > //j + 3i/4, the Bethe ansatz equations f|70[) and (f?2")) can be written as 



Aj + i/4 



n iV M 2 



Aj - i/4 



Aj — /ii + i/2 



JJ A, - Mi - i/2 11 A, - A, - i/2 



Mi 

n 



Aj - Aj + i/2 



Mi . . , M-2 . 

yr fXj - Aj + z/2 _ _ tt H - fM + I 

r=i ^ ~ A « ~ i / 2 i=i & ~ Mi - * ' 

From the relation (|7Tj) . we obtain the eigenvalue i? of the Hamiltonian H {T]) as 

Mi 



^_ 9J/dlnT(A) 



9A 



11 



A=0 



9Jtt 



99 



(73) 



(74) 



j=i 



where a t (x) = t/[Tr(x 2 + i 2 )] and Aj should satisfy the Bethe ansatz equations (|73[) . The momentum P is pi. [2r| 

(75) 



P = fcj mod 27r, 

3=1 

where kj are parameterized by e zkj = (Aj — «/4)/(Aj + i/4). 

The basis of Hilbert space Vi are |3/2)., |l/2)-, — |l/2)- and — |3/2).. Thus the vacuum state is a ferromagnetic 
state |0) = 1 3/2) 1 <8>|3/2) 2 <8>- ■ -<8>\3/2) N . The M1+M2 is the total number of flipped spins from this ferromagnetic 
state, so the total spin along the z-component and magnetization are 

S z = 3N/2- (M 1 +M 2 ), (76) 
m = 3/2 - (Mi + M 2 )/N = 3/2 - n x - n^. (77) 

From Eq. (|74p . if J < 0, the A will lead to addition to the energy. Therefore, the ground state of this case is 
the ferromagnetic state |0). While if J > 0, the A will contribute a negative value to the energy, so the ground 
state might be a anti-ferromagnetic state. In order to obtain the ground state configuration of the system, we 
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IV. THERMODYNAMICS 



Now, we solve the Bethe ansatz equations (|73)) . In the thermodynamic limit, both the total particle number 
and the system size tend to infinity while the ratio N/L keeps a non-zero constant. The Bethe ansatz equations 
can have the complex solutions, i.e., strings. Checking of the Bethe ansatz equations (|73p in detail, we find that 
the string hypothesis of the Bethe ansatz equations (|75|) is 

A« = A£ + |(n + l-2j), j = l,...,n, 

AC = AC + \{n + 1 - 2.?), j = 1, ■ ■ • , m, (78) 

where XI and /z™ are the real parts of the n-string of A and the m-string of /i, respectively. From Eq. (fTTj) . we 
see that the contributions of n-string of A to the energy is 

e«(A£ l ) = -\j™ n/ ,{\ n z ). (79) 

The energy E and the momentum P are 

E = J2e n (\ n z )-Y JN = -l J *Y, a n/*W)-Y JN > (80) 

n,z n,Z 

p = E 7 ' mod 27r > ( 81 ) 

n,z 

where [x — it)/(x + it) = —e l6t ^ and 9t(x) = 2 arctan(x/i). 

Substituting the string solutions into the Bethe ansatz equations (|73[) . we have 



e -*mn {K ) jrj e -iB^„(A»-^) = (_)<5x- JJ e -^ 4 m ,n(A?-A™) ; (82) 
J] e - iB «.™^- x ™) - (-) 5 «" 1] (83) 

where 



B m ,ni x ) — ^i(2m+n-l) + ^|(2m+n-3) + ^i(|2m-n|+l) ' ( 84 ) 
•Kn,n{ x ) = ^i(m+n) + 2 $i(m+n-2) "I + 2 ^i(|m-n|) +^i(| TO -n|)) (85) 



and <5 gives the phase shifts n or according to the numbers of production terms. Taking the logarithm of (1831) . 
we have 

2nl* z = N9n (A") + E B^,,«(A; - - E ( A " " 

where 7^ 2 and 7^ are the integers or half-integers which determine the eigenstates. A set of {7^ z , 7^ } satisfying 
(|86|) gives a highest weight eigenstate of Hamiltonian |T]). The momentum (|75|) can be written in terms of 7„ 2 
and I? ny 

P = 2tt1^(7 A; +7 a? ) mod27r. (87) 

In the thermodynamic limit, the summations become integrations. Denoting rj n and a m as the densities of A 
n-strings and /j, m-strings in the thermodynamic limit, and ?7„ and cr^ as the corresponding densities of holes, 
the Bethe ansatz equations read 



N 



L = e n/4 (X) + £ Bt, n * <r m (A) - E * */ m (A), 

m m 

= E^, m *™-E^,»*^ m (M), 
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where the densities of A™ and /x" are denoted as r] n (X z ) and r] n (fj, z ), respectively, Is are the quantum numbers, 
and the operator * is defined by 



a t (x-y)f(y)dy=[t]*f(x). (89) 
Taking the differentials of Eq. we obtain the integral form of the Bethe ansatz equations 

r,l{k) = a n/4 (k) + £ Al„ * * m (k) Yl A kn * V m (k), 

ni—1 ni—1 

aUk) = J2 Km * V m (k) - £ * a m (k). (90) 



m— 1 



where ^„ and „ arc the integral operators, A l m n = a( m+ „)/ t +2a( m+ „_ 2 )/ t H h2a(| TO _ n |_ 2 )/t+a(| m _„|)/ t , 

n = a (2m+ii-i)/(+ a (2m+n-3)/(+' ' ' + a ( \2m-n\+i)/t ■ In the derivation, we have used the relation (d/dx)6 t (x) = 
2ira,t(x). 

At temperature T, the Gibbs free energy of the system (1) with an external magnetic field h reads 

F = E-h(%N -M-l-Mv) -TS, (91) 



2 



where 



My = N £ n J r, n (A)dA, M 2 = N^m J a m (/z)d/z (92) 
£ = --ttJJV^ / dAa„ /4 (A)rf (A) - -JV, (93) 

n 

S = NJ2 [ dA[(r?„ + ^) ln(r?„ + t£) - 7? n In r?„ - 77* In t£] 

+N J2 M(Pm + a ra) M^n + o£) ~ <7 m hi <7 m - O h m In ] . (94) 
m •* 

Minimizing the Gibbs free energy at the thermal equilibrium, we obtain the following thermodynamic Bethe 
ansatz equations 

Inm = —ttJ^- + G 4 * ln[(l + + p" 1 )" 1 ], 



In ^ = G 2 * In --it/?, — u ~ 7^ * ln (! + % 1 )> 



In ?]neeven = G 4 * In — , (95) 

1 + <T n/2 

ln7)„ eocW = G 4 * ln[(l + ry„_i)(l + f/ n+1 )}, 

i ~ /-< i (1 + ^m-l)(l + ^m+l) G*2 . , — h 
lna m = G 2 * In -i — — ^- \ - — * ln(l + 



lnr?„ /i ln<r m ft, 
lim = — , urn = — , 

n — >oo 77 T m — >oc 777 T 

where fj n — Vn/Vn- &m ~ ^m/^mi G„ = a 1 /„/(a + a 2 /„) and qq = S(x). In the derivation, we have used the 
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following relations 



[n + m\ 

A l,m 



[n] [to] 



.4 



n+l,m 
t 



) = 5 m ,n, [n > 1), 



m,n 
t 

m,l 
t 

n,m 
t 



G t {B™ n - 1 + B^ n+1 ) = S 2m , n , (n > 1), 



G t B?' 2 = 0, 

n— l,m 



G t / 2 (-B. 



n+l,m 



) — (<^2m-l,n + ^2m+l,ra)G t /2 

+ 8 2m , n G t / 2 /G tl (n > 1), 



5 t 1,m - G t/2 5 t 2 ' m = (J 1)n + 62, n )G t/a (n > 1) + 5 2 , n G t/2 /G t , 



II 


+ 


r 


, rn 


"71" 




t 




.7. 


n 


+ 


r 


D 771,71 


"71" 




t 




.1. 


n 


+ 


i" 




n" 


t 


.1. 



.4 



n+l,m 



m,n+l 



n+l,m 



0, (m < n) 

-a m / t /Gt, (to > n + 1) ' 
0, (2m < n) 

-a 2m /t/Gt, (2m > n + 1) 
0, (to < 2n) 

-a(m+i)/t! ( m = 2n+ 1) 
-a m / t G t , (m > 2n + 2) 



(96) 



V. GROUND STATE 



In order to obtain the ground state of the system (1), we first define the dressed energies of n-string A and 
m-string (i 



CnW = T In fj n (\), ? m (/i) = Thi(7„(/j). 
From Eqs. ([50)) and (JHSJ), we find that the dressed energies should satisfy 

Cn(A) = -§7rJa„ /4 (A) + ftn - £ * M 1 + e- ? '"< A >/ T ] 



(97) 



*,(*:) = ftn - ^ B^ TO * ln[l + e"«" 



(A)/Ti 



m=l 

Ei^„*ln[l + e--W/ T ]. 



(98) 
(99) 



m—1 



In the thermodynamic limit, we obtain the following thermodynamic Bethe ansatz equations for the dressed 
energy 

Ci = -^nJG 4 (\) + TGi * ln[l + e^' T ], 

ft = TG 2 * {ln[l + e^/ T ] - ln[l + e -^ T ] - ln[l + e^^]} 

,G 2 



_T^*ln(l + e^( A )/ T ), 

Cnecvcn = Td * {hi[l + e<-^ T + ln[l + e^' T ] - ln[l + e"^/^ 

C„ eo dd = TG 4 * {ln[l + e^' T + ln[l + e^' T ]}, 

W = TG 2 * [ln(l + e^- l/T ) + ln(l + e ?m + l/T ) - ln(l + e -6»+i/T) 

- ln(l + e-C--i/ T )] - * ln(l + e -^™/ T ), 

G4 

,. Cn , ,. ?m , 

lim — = /i, lim — = ft. 

n— >oc 72 m^oo TO 



(100) 



The dressed energies can be divided into two parts e + (fc) and e (A;) 



e(fc) if e(fc) < 0, 
if e(fc) > 0; " 



(101) 



The ground state string distribution of the system can be obtained by taking the limit of T — > and h — > 0. 

a 17 ^ 1 in, WWTi 1 r-i , ,+ flWT rlll/Ti 1 i r, , -tflVTl i r, , 
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e -e+(A)/r e -«-(A)/T]_ Whcn T ^ Q) both e -e+(A)/T and e £-(A)/T tend tQ one Xhug m [l +e e(A)/T] = m [l + e £ + (A)/T] 

and ln[l + e~ e( - x V T } = ln[l + e - £ "( A )/ T ]. 

If J < 0, the ground state of the system is the ferromagnetic state |0) = ®jLi 1 3 / 2} . It is easy to understand 
from the eigen energy (|80[) for a n-string A give a positive contribution e n {\ z ) to the eigen energy. In the 
ground state, the total spin along z-direction and the magnetization are S z = 3N/2 and m = 3/2, respectively. 
It is a highest weight representation of the Yang-algebra The ground state energy and momentum are 
E = -99JN/8 and P = 0, respectively. 

15 
10 

5 



-J 

Uj -5 
-10 
-15 
-20 
-25 




FIG. 1: The ground state energy of the system. There is a quantum phase transition at the critical point J = 0. 
system is in the ferromagnetic phase if J < and is in the antiferromagnetic phase if J > 0. 



The 



If J > 0, wc find that some As are real while the others form 2-strings, and the /js are real. Such a ground 
state configuration is quite different from that of the SU(A) Sutherland model where there is no string or spin 
bound state in the ground state. In the present SO(5) case, part of the spectral parameters form 2-strings which 
heavily affect the spin excitations as we shall show below. To make the ground state energy lowest, all these 
strings arc filled up and no holes left. This can be understood from the entropy S of the system. The completely 
filled string configurations contribute zero entropy, and make the system in a most stable sate. 

From Eqs. ([501) . the densities 7$ and cr™ satisfy the following integral equations 



Po(A) = <?(A) +K *p (A), 
where p (X) = [r/oi (A), 7702(A), 001(A)]*, g(X) = [a 1/4 (A), a 2 / 4 (A), 0]' and 



K = 



-<Xi 

2 

ai 
a i 



-S3 — ai 
4 4 

— ai — 2ai 
a,3 + ai 




Taking the inverse of Eq. (|102[) . we obtain 

po (A) =F*g{X), 
where F = 1/(1- K). The solution of (fTM)) is 



(102) 



(103) 



(104) 



7701(A) = cosIi(ttA), 
1 

187r 

001(A) = cosh(7rA/3)/3. 



%2(A) 



-csch(2-7rA) 



47r\/3 sinh 



4ttA 



12ttA 



(105) 



Because all the density functions 



are 



even, the ground state string configuration 



Uo 2 (M 2) )L UoVl^)} is symmetric around the origin. Taking the integration, we obtain the 



densities of strings = M^/N, 



7} 



Ml/N and nfc = M\jN as 



'01 



'02 



(106) 



After some derivations, we find M\ = N and M2 = N/2, which mean that the total spin of the ground state is 
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FIG. 2: The densities of »? 1 (A), rj (A) and o" 1 (A) of the ground state. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 

V 

FIG. 3: The dressed energies C 1 (A), C 2 (A) and £ (A) of the ground state with J = 2/9ir. 
and the momentum are 

9 99 

= --wJNF(0) - — JN, F{\) =a 1/4 * 7701(A) + a 1/2 *7fo 2 (A), (107) 
P A = 27T-i^(/ A y+7 A ») mod 2^ = 0. (108) 

71,2 

From Eq. (J95J), the ground state dressed energies satisfy the following equations 

e(X) = -^nJF * .9(A) = -^rJpb(A), (109) 

where e(A) = [Ci(A), <^i(A)] T . Recalling the definition of the dressed energy e = T\n(p h /p), we see that the 

dressed energies are negative in the limit of T — > 0, which means that the corresponding strings are completely 
filled. 



VI. ELEMENTARY EXCITATIONS 

Based on the ground state configuration, the elementary excitations of the system can be studied exactly. 
For the ferromagnetic case (J < 0), elementary excitations are spin waves with dispersion relation AE = 
— 18( J/n) cos 2 (nAP/2). The excitations in anti-ferromagnetic sector (J > 0) are somehow complicated. In the 
language of Bethe ansatz, these can be described by the changes of the string distributions. These excitations 
are very different from those of the SU(A) model, for the ground state configuration contains 2-strings. In fact, 
the spin excitations can be described by adding some holes or high strings into the ground state configuration. 

The holes and extra strings lead to redistributions of As and [is. Formally, the extra strings contribute nothing 
to the energy because the contribution of such strings is exactly canceled by the rearrangement of the ground 
state distribution though they do contribute to the spin quanta carried by the excitations. 
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FIG. 4: The single-hole excitations of the system. Here J = 2/9n, AE and P are the energy and the momentum carried 
by a single hole, respectively. The dotted dashed line is the single-hole excitation of real A. The solid line is that of A 
2-string and the dashed line is that of real p. 



The excited states can be determined by the integral Bethe ansatz equations with holes and high strings in 
the A , A 2 and p 1 sectors 

p(X)=F*tf(X)-ff l (X)}, (110) 

where p(X) = [771(A), 772(A), cri (A)] T and /^(A) = [77^ (A), rfe (A), a\ (A)] T . In the thermodynamic limit the density 
of holes are 

ma -i 
1=1 

where m 1: 2,3 represent the numbers of holes in real A sea, in 2-string A sea and in real p sea, i/*, and i/| are the 
positions of the corresponding charges and holes. The excitations lead to the redistributions of densities po(X), 

p(X) + f? l (X) - po(k) = K * [/7(A) - p (k)}. (112) 

Denoting the charges of the densities as Ap(k) = p(X) — po(X), which satisfies 

Ap(A) = -F*p h (X). (113) 

Form Eq. ([93]) . the excited energy is 

3 m a 

AE = E-E Q ^ = -Nj2^*p h a (X)\ x=a = -J2J2 e ^<y ( 114 ) 

a a—1 i—1 

Thus the excited energies are the summation of dressed energies carried by the holes with a inverse sign. The 
excited momentum is 

3 m a „ v * 

AP = P - P A = / " A}a(A)dA mod 2tt. (115) 

a=l i=a J° 

The spin quanta carried by the spin excitation is 

S = |m 2 + 2m 3 + ^)(2 - l)m xm + ^(1 - t)m^ t) , (116) 

Z>3 t>2 

where m A (i) and m^t) are the numbers of A /-strings and the p i-strings formed in the excitations, respectively. 

Because the energies, momenta and spins of the holes are additive, the thermodynamic behaviors of the system 
are mainly determined by the dispersion relations of the individual holes, which are shown in Fig[3] From the 
FigJ21 we find that the single A 2-string hole carries the lowest energy with spin 3/2, which is named as heavy 
spinon here. These heavy spinons dominate the low temperature thermodynamics of the system. Surprisingly, 
the holes in the real A sector carry zero spin, corresponding to a new kind of neutral spin excitations. The spin 
quanta carried by each p hole are 2. Different from the SU(A) model, the A 2-string heavy spinons cover one 

pji n •n- ■ /_.. t^' m ml.' . . * . _l _ j l l.i. .j. 11 ii 1. ...ii. ' . .11 l. 
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TABLE I: Some possible configurations of the excitations. Here, AM a = AMI an d m? are basic excitations. m A ; and 
mj 1 are the high string excitations. Some possible elementary excitations mi = mi+nij , m.2 = mjj+m^ , 1113 = mi + 2m3 , 
ni4 = mi + mg + 2mj and ms = mi + mg + mj are shown. Here a x b means the number of ^-strings added is a. 



m 


mi rri2 mz m x (i) m^t) 


AA/1 1 AMi AAf! 


A Mi AM 2 


AS Z 


m? 
mi 
mi 
mi 
mi 
mi 


2 
10 10 
1 2 
2 10 
4 
2 


-2 1 
-1 -1 

-1 -1 

1 -2 -2 

2 -3 -2 
-1 -2 




-1 -1 
-2 -1 
-3 -2 
-4 -2 
-2 -2 



2 
3 
5 
6 
4 


m^ 
m? 


1 x I 
lxt 


0-10 
0-1 


2 - 1 
t-1 


-I 
-t 


mi 
m 2 
m 3 
m.4 
m 5 


10 1 1x2 
1 2 1x4 
1 2 2x3 
4 1x62x2 
4 1x61x3 


-1 -2 
-2 -1 
-3 -1 
2 -4 -4 
2 -4 -3 


-1 
-1 
-1 




1 
1 
1 





The numbers of holes and strings added are not independent but satisfy some constraints determined by the 
Bcthe ansatz equations (|90j) . 

AMI = -mi + ^m 2 , 

AMI = -mi - -m 2 - -m 3 - ^ m x m , (117) 

l>3 

AMI = ~^ m 2 - TO 3 ~ J! «Vt) , 

t>2 

where 7711,2,3 are the numbers of holes in the real A, 2-string A and real /x-sea, m^m and rriM) are the number of 
Z-string in the rapidity A and t-string in the rapidity ix, respectively. Thus AM[ 2 are some integers which indicate 
the number changes of A, /1 Z-strings. For convenience, we define AM a = J^i AM* and denote these excitations as 
m = [(mi, m2, 7773), (m^, m l x , • • ■ ), (m^, m^ ,•••)]- As we mentioned above, the excited momenta, excited energies 
and the spins are additive, AP m+m / = AP m + AP m , , AE m+m , = AE m + AE m , and S^ +m , = + S z m , . The 
number changes of AM% are also additive 

AM„ fc m+m , = AM 6 m + AM a b m , , AA/ am+m , = AM am + AM am , . (118) 

Some possible hole configurations are listed in Table Q] We denote [(a, b, c), (0, 0, • • • ), (0, 0, • • • )] = [(a, b, c)] 
for short. We find that the excitations m° = [(2,0,0)], = [(1,0,1)], m° = [(1,2,0)], = [(0,0,2)], 
mi = [(0,2,1)] and mi = [(0,4,0)] are the basic excitations. The additional single high strings m A ; are not 
independent. They can form the excitations with m°, such as mi, m.2, m.3, m.4 and ms. 

The low-lying excitations m° are shown in Fig[2] The simplest spin excitation is a real A hole-pair m° that 
is two-neutral spinon excitation, corresponding to the two domain walls of a single excited domain (Fig. EJa)). 
The A 2-string hole pair can not exist independently. They must be associated with a neutral spinon, i. e. m§ 
(Fig. EJc)). Accompanied by a real /j hole, the A 2-string hole pair is also a possible excitation m" (Fig. E] (e)). 
Further, if we add a A 4-string into a 2-string hole pair and a real A hole in the A-sea, the total spin of this 
excitation is 1 (rri2 in Table. []}. ln this case each of the A 2-string holes carries a spin 1/2. Such a dressed hole is 
quite similar to the ordinary spinon, and it is named as dressed spinon here. Four A 2-string holes mi may exist 
independently (Fig. [5]). If we put further one A 6-string and one fi 3-strings into this four hole configuration 
(1115 in Table. we obtain the 50(5) spin singlet excitation. The simplest excitation in the fi sector is a pair 
of real /i holes mi (Fig. EJd)). This excitation is quite similar to a real A hole pair but each of the real /x hole 
carries a spin 2. Joint pair of a real A hole and a real [i hole m§ may also happen as shown in Table. |TJ Other 
kinds of spin excitations such as mi, 1113 and m 4 can be constructed similarly. 
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-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



A/Vji A/Vji A/Vji 



FIG. 5: The low-lying excitations of the system. Here J = 2tv/9, AE and AP are the energy and the momentum carried 
by the excitation, respectively, (a) m? = [(200)]; (b) m§ = [(101)]; (c) = [(120)]; (d) = [(002)]; (e) m 1 ] = [(021)]; 
(f) mg = [(040)]. 



VII. CONCLUSION 

In conclusion, we propose an integrable spin-3/2 chain model with 5*0(5) symmetry. By using the nested 
quantum inverse scattering method, we obtain the exact solutions of the system. Different from the 57/(4) 
integrable spin chain, there only exist three conserved quantities. Based on the exact solutions, the ground state 
and thermodynamic properties of the system are discussed. Several new kinds of spin excitations such as the 
neutral spin excitations, heavy spinons and dressed spinons are found. 
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